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Abstract 

We derive the model equations describing the thermoacoustic res- 
onator, that is, an acoustical resonator containing a viscous medium in- 
side. Previous studies on this system have addressed this sytem in the 
frame of the plane-wave approximation, we extend the previous model to 
by considering spatial effects in a large aperture resonator. This model 
exhibits pattern formation and localized structures scenario. 

1 Derivation of the model 

The physical system we consider is an acoustic resonator composed of two par- 
allel walls containing a viscous fluid medium inside (e.g. glycerine). The case 
of sound beam propagating in viscous media is characterized by a strong ab- 
sorption, being the sound velocity dependent on the fluid temperature, which 
results in an addicional nonlinearity mechanism from thermal origin. 

1.1 Sound propagation in a viscous medium 

The propagation of sound in a viscous, heat-conducting medium is described by 
the equations of mechanics of continuum media [T], i.e. the continuity equation 

§^ + PoV-v = 0, (1) 
the momentum transfer equation 

^ + jVp = p^iiV\ + (^fiB + ^) V (V • v) , (2) 
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and the heat transport equation 



dT dv^ ( dv^ dvk 2 dvi\ 2 

where p and p represent, respectively, the density and acoustic pressure devia- 
tions from the equihbrium state (denoted by pg and po), v is the fluid particle 
velocity. The constants /j,, /i^)'* ^.nd Cp are the coefficients of shear and bulk 
viscosity, the thermal conductivity and the specific heat at constant pressure 
of the fluid, respectively. In Eq.Q, T denotes the perturbation of the medium 
temperature due to the acoustic wave, and Xi the spatial coordinates. 

Equations (P)-© must be complemented with the equation of state relat- 
ing the medium density with the other variables. For a linear medium (small 
amplitudes) it obeys the simple form p — Cgp. In a highly absorbing medium, 
the thermal effects are dominant and the hydrodynamic nonlinearity can be 
neglected. Furthermore, the speed of sound depends on the average value (over 
one period) of the temperature Tav = Tq + T, where Tq is the temperature at 
equilibrium, so the equation of state takes the form p = p/cq{T). Expanding at 
leading order it follows 

P-ip-^pT- (4) 

where 7^ = (2/co) {m)T=T,, 

We assume that the fields can be described by quasi-plane waves, i.e. they 
take the form 

p = ip(r,i)e*-* + c.c., (5a) 
V = iu(r,i)e*'^*-f c.c, (5b) 
T = ir(r,i)e*'"* + c.c. + r(r,t), (5c) 

where a slowly changing temperature T, responsible for the thermal self-action 
of the waves on w = Vz, is added. Taking into account the slowly varying 
envelope condition 

dX, , dX, 

^«kX.,^«.X., (6) 

and the linear impedance relation relation v — p/pqCq, the equation for T is 
obtained by averaging of Eq. ([3]) 



dt ^ 2plcl ' ' ' 

where T represents the average temperature perturbation caused by the wave 
energy dissipation. In Eq.© the new parameter b is defined as b — (p^ + 
4/x/3)/pqCp, and x = i^/Po'^p the thermal diffusivity. 
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On the other side, substituting Eqs.® in Eqs.lH)) and ^ it follows that, in 
the parabolic approximation, the pressure field is described by \1\ 

dP dP ic?. , 7n - 

coj- + ^ + Tr^lP + coaP = tuj^PT (8 
dz ot 2uj 2 

where a = buj"^ /2pqCq and is a laplacian operator acting on transverse coor- 
dinates. 

1.2 Introducing the resonator 
1.2.1 Modal equations 

In a resonator, the reflection at the boundaries imply the presence of counter- 
propagating traveling waves, and the acoustic field p can be expressed in terms 
of forward and backward components 

P = PF(r,i)e*('^*~''^' +Ps(r,t)e^("'*+'=^'. (9) 

whose amplitudes are assumed to vary slowly in space and time (compared with 
the evolution described by the exponential factors). 

Substitution of expansion ([9]) into Eq. ^ , leads to the following temperature 
distribution 

T = T^(r, t) + Tgir, i)e*2/c. ^ f)e-'2fc.^ (^q) 

where Tq is the homogeneous component and Ti is the amplitudes of the thermal 
grating induced by the acoustic wave. The acoustic resonator containing a 
viscous fluid medium inside (e.g. glycerine), has been previously addressed [2]in 
the frame of plane-wave approximation, we extend now the previous model by 
considering the sound diffraction and temerature diffusion. 

Projecting on the different longitudinal modes, the following equations for 
the amplitudes are obtained: 

iuj^ {ThPf + T^Pb) , (11a) 
^u;^(ThPB+TgPF), (lib) 



dPp 


dPp 




dt 


dz 




dPs 


dPs 




dt 


dz 


2w 



CQaPf 



coqPb = 



^-xVir. = ^(\Pf\' + \Pb\'), (12a) 



dt - - " 2plc, 

dT„ ... Pb 

a 



^+4xkX~X^lT, = ^PfPb, (12b) 



where a = w {fig + 4^/3) / {2pqCq). 
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1.2.2 Boundary Conditions. 

We consider an acoustic resonator with length L, filled with a viscous medium, 
and bounded by reflecting surfaces, with intensity reflection coefficients TZ, lo- 
cated at z = and z = L. We assume that one of the surfaces, vibrating 
at frequency uj, acts as ultrasonic source. In accordance with ([H]), boundary 
conditions satisfy 

PF(r±,z = 0;t) = V^PB(ri,2 = 0;i) + P™, (13) 
Pb{v±,z ^ L,t)^VnPF{r±,z^L;t)e-'\ (14) 

where the detuning parameter S has been defined as 

S = 2m7: - 2kL = 21 ^" ~ ^ , (15) 
Co 

where m is an integer, selected in order to be lOc the cavity frequency that lies 
nearest to the driving frequency uj. 

1.2.3 Mean Field Limit. 

When the reflectivity of the boundary walls is high, TZ ^ 1, the resulting model 
can be greatly symplified adopting the mean field limit. Under this assumption, 
the two counterpropagating fields are approximately constant and equal along 
the cavity (i.e. there is no spatial modulation of the field envelopes), and one 
can assume that 

Pf(r±,i)=PB(rj.,i) = P(r^,i) (16) 
where the overbar means the spatial average over the longitudinal coordinate z, 

L 

I^\jdzf{z). (17) 



We also consider that the frequency detuning between the driving and a 
cavity mode is small, 2L (ujc — w) /c << tt, which allows to write the boundary 
condition Eq.((T4]) as 

PB{ruL,t) « v^PF(r±,L,t)(l-i,5). (18) 

Averaging Eqs. lfTTjl . with 

i hPf^b,-.- Pf^^ iL)-Pf^, (0) 
lJ dz L 



iLj^ (ThPf + T^Pb), 
tiu^ {ThPb + TgPp) , 



leads to 

dPp 

BPb 



f [Pf {L) 



^F(0)]-^^ViP^. 



--|[P.(X)-P.(0)]-.|-ViP. 



caapB = 
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they can be reduced just to one equation 



(19) 

where the averaged temperature components as defined as Ti = CTi and we have 
used the property {fg) = / 5 [3], valid in the mean field case. The transmitted 
pressure can be assumed to be 

PT(t) = VTP(r±;t) (20) 

and Eg. p^ is written as 

dtP = -r(l + t6)P + P„ + »^VlP + tuj^ Tn + P, (21) 



where the parameters 

V = n. r.n A- 1 

2L 



aco+T^, (22) 



(23) 



r 

2L 

have been introduced. As well, under the same assumptions, Eqs. p2)) read 



P^n = ^Prn, (24) 



xVln - — \Pr ; (25a) 
^ + Akhf, - xVif, ^ IP^, (25b) 



1.3 The original model 

The mean-field model for the thermoacoustic resonator can be written as 

2 

9tP = ~^{T + 26)P+^P,^ + t^VlP~coaP + t^{n+fgl^^) 



2L' ' 2L 2uj ^ 2 

k^b 



dt% = _4fc2;^r, + xVir, + -^|P|', (26c) 

where Tg has already been assumed to be real. We note that the model is 
highly reminiscent of that for a coherently driven optical cavity filled with a 
non- instantaneous (and nonlocal) Kerr medium [4]. Unlike its optical analog, in 
this thermoacoustic model there are two nonlinearities (owed to Th and to Tg); 
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notice however that for 7/j = 4fc^x both nonhnearities behave as a single one (a 
single temperature field Th + Tg can be defined so that its evolution and that 
of the pressure just depend on themselves in that case). Defining the following 
quantities 

c T 

tp^ = -jj^ + coa, A = tp {ujc -uj), Tp = 7htp, Tg = (27) 



V 47hPoCo 
the model equations, Eqs. ((26| become 




2V^p^ ^ ^ jj^j.^^ C ^ (29) 



TpdrP = -{l + iA)P + Pin + iy^P + i{H + G)P, (30a) 
= -H + DV'^H + 2\Pf , (30b) 
a^G = -T-^G + DW^G+\P\\ (30c) 

We consider a resonator with high quality plates (T ~ 0.1), separated by L = 5 
cm, driven at a frequency / = 2MHz and containing glycerine at 10°C. Un- 
der these conditions the medium parameters are cq = 2 x lO^m s^^,q;o = 10 
m-i, po = 1-2 X lO^kg m-^, Cp = 4 x lO^J kg-^K~^, a = 10"^ ^ ^ q.S 

W m-iK-i(x = lO-^m^s-i). In this case tp = 2 x 10"^ s, = 6 x 10"^ s, 
and our length unit is = 2 mm. For a resonator with a large Fresnel number, 
the relaxation of the homogeneous the temperature is mainly due to the heat 
flux through the boundaries, and can be estimated from the Newton's colling 
law as tfi ~ 10^ s. Then, under usual conditions we get the normalized decay 
times Tp ^ 10~^, and Tg ~ 10~^, and the diffusion constant D ^ lO'^.Note that 
for Tg = 1 a single temperature field {H + G) can be defined, whose evolution 
depends just on P and vice versa. This model exhibits a rich spatiotemporal 
dynamics and exhibits pattern formation and localized structure formed in the 
transverse cross-section of the cavity. Some of these results are studied in [5] . 
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